Chapter 35 



Quantum scattering 



Scattering is easier than gathering. 
— Irish proverb 

(A. Wirzba, P. Cvitanovic and N. Whelan) 

So FAR the trace formulas have been derived assuming that the system under 
consideration is bound. As we shall now see, we are in luck - the semiclas- 
sics of bound systems is all we need to understand the semiclassics for open, 
scattering systems as well. We start by a brief review of the quantum theory of 
elastic scattering of a point particle from a (repulsive) potential, and then develop 
the connection to the standard Gutzwiller theory for bound systems. We do this 
in two steps - first, a heuristic derivation which helps us understand in what sense 
density of states is "density," and then we sketch a general derivation of the cen- 
tral result of the spectral theory of quantum scattering, the Krein-Friedel-Lloyd 
formula. The end result is that we establish a connection between the scatter- 
ing resonances (both positions and widths) of an open quantum system and the 
poles of the trace of the Green's function, which we learned to analyze in earlier 
chapters. 



35.1 Density of states 

For a scattering problem the density of states (31.16) appear- ill defined since for- 
mulas such as (34.6) involve integration over infinite spatial extent. What we will 
now show is that a quantity that makes sense physically is the difference of two 
densities - the first with the scatterer present and the second with the scatterer 
absent. 

In non-relativistic dynamics the relative motion can be separated from the 
center-of-mass motion. Therefore the elastic scattering of two particles can be 
treated as the scattering of one particle from a static potential V{q). We will study 
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the scattering of a point-particle of (reduced) mass m by a short-range potential 
V{q), excluding inter alia the Coulomb potential. (The Coulomb potential decays 
slowly as a function of q so that various asymptotic approximations which apply 
to general potentials fail for it.) Although we can choose the spatial coordinate 
frame freely, it is advisable to place its origin somewhere near the geometrical 
center of the potential. The scattering problem is solved, if a scattering solution 
to the time-independent Schrodinger equation (??) 



can be constructed. Here E is the energy, p - fik the initial momentum of the 
particle, and k the corresponding wave vector. 

When the argument r - \q\ of the wave function is large compared to the typ- 
ical size a of the scattering region, the Schrodinger equation effectively becomes 
a free particle equation because of the short-range nature of the potential. In the 
asymptotic domain r » a, the solution (fi^iq) of (35.1) can be written as superpo- 
sition of ingoing and outgoing solutions of the free particle Schrodinger equation 
for fixed angular momentum: 

(j){q) = A(l)'-~\q) + B(p^^^{q) , (-1- boundary conditions) , 

where in 1-dimensional problems (p'-~\q), 4>^'^\q) are the "left," "right" moving 
plane waves, and in higher-dimensional scattering problems the "incoming," "out- 
going" radial waves, with the constant matrices A, B fixed by the boundary con- 
ditions. What are the boundary conditions? The scatterer can modify only the 
outgoing waves (see figure 35.1), since the incoming ones, by definition, have yet 
to encounter the scattering region. This defines the quantum mechanical scattering 
matrix, or the S matrix 



All scattering effects ai^e incorporated in the deviation of S from the unit matrix, 
the transition matrix T 



For concreteness, we have specialized to two dimensions, although the final for- 
mula is true for arbitrary dimensions. The indices m and m' are the angular mo- 
menta quantum numbers for the incoming and outgoing state of the scattering 
wave function, labeling the 5 -matrix elements S,nm'- More generally, given a set 
of quantum numbers yS, y, the S matrix is a collection S py of transition amplitudes 
13 ^ y normalized such that 15^3^1^ is the probability of the {3 ^ y transition. The 




(35.1) 



Mr) = 4>m V) + S mm'0^,V) . 



(35.2) 



S - 1 - /T . 



(35.3) 



scattering - 29dec2004 



ChaosBook.org versionI4.4.1, Apr 21 2013 



CHAPTER 35. QUANTUM SCATTERING 



673 



Figure 35.1: (a) Incoming spherical waves run- 
ning into an obstacle, (b) Superposition of outgo- 
ing spherical waves scattered from an obstacle. 



(a) 



(b) 



total probability that the ingoing state /3 ends up in some outgoing state must add 
up to unity 



(35.4) 



so the S matrix is unitary: S^S = SS^ = 1. 

We have already encountered a solution to the 2-dimensional problem; free 
particle propagation Green's function (33.48) is a radial solution, given in terms 
of the Hankel function 



Go{r,0,E) = -^Hl^'-\kr), 



The mth angular mo- 



m \Q) Hl^\kr), and given a potential 



where we have used S Q(r, 0, E) /H = kr for the action 
mentum eigenfunction is proportional to 
V{q) we can in principle compute the infinity of matrix elements 5 We will 
not need much information about H^^{kr), other than that for large r its asymptotic 
form is 



oc e 



tikr 



In general, the potential V(q) is not radially symmetric and (35.1) has to be 
solved numerically, by explicit integration, or by diagonalizing a large matrix in 
a specific basis. To simplify things a bit, we assume for the time being that a 
radially symmetric scatterer is centered at the origin; the final formula will be 
true for arbitrary asymmetric potentials. Then the solutions of the Schiodinger 
equation (??) are separable, (pmiq) = (p{r)e'"^^ , r = \q\, the scattering matrix cannot 
mix different angular momentum eigenstates, and S is diagonal in the radial basis 
(35.2) with matrix elements given by 



S,n{k) 



(35.5) 
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The matrix is unitary so in a diagonal basis all entries are pure phases. This means 
that an incoming state of the form H\^\kr)e^'"^ gets scattered into an outgoing state 
of the form S ,„{k)H^^\kr)e^'"^ , where hI^\z) are incoming and outgoing Hankel 
functions respectively. We now embed the scatterer in a infinite cylindrical well 
of radius R, and will later take R ^ oo. Angular momentum is still conserved so 
that each eigenstate of this (now bound) problem corresponds to some value of m. 
For large r :s> a each eigenstate is of the asymptotically free form 



where • • • is a common prefactor, and;^„, = mnll+njA is an annoying phase factor 
from the asymptotic expansion of the Hankel functions that will play no role in 
what follows. 

The state (35.6) must satisfy the external boundary condition that it vanish at 
r = R. This implies the quantization condition 

knR + 5,n{kn) - Xm = TT (ft + 12) . 

We now ask for the difference in the eigenvalues of two consecutive states of 
fixed m. Since R is large, the density of states is high, and the phase 6,nik) does 
not change much over such a small interval. Therefore, to leading order we can 
include the effect of the change of the phase on state « + 1 by Taylor expanding, is 

kn+lR + 5r„{kn) + (^„+l - kn)5'^{kn) - Xm ~ TT + n(n + 12) . 

Taking the difference of the two equations we obtain Ak k n{R + 6'„{k)yK This 
is the eigenvalue spacing which we now interpret as the inverse of the density of 
states within m angular momentum sbuspace 



The R term is essentially the \ - d Weyl term (34.8), appropriate to I - d radial 
quantization. For large R, the dominant behavior is given by the size of the circular 
enclosure with a correction in terms of the derivative of the scattering phase shift, 
approximation accurate to order l/R. However, not all is well: the area under 
consideration tends to infinity. We regularize this by subtracting from the result 
from the free particle density of states do{k), for the same size container, but this 
time without any scatterer, figure 35.2. We also sum over all m values so that 



Mr) ^ e'"''>[s„^{k)H^^\kr)+H^-\kr)) 



X ■ ■ ■ cos{kr + 6,„{k) - Xm) , 



(35.6) 



dm{k) x-{R + 6'„m . 




m 



m 




(35.7) 
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Figure 35.2: The "difference" of two bounded refer- •, b 

ence systems, one with and one without the scattering 

system. ^ ^ 



The first line follows from the definition of the phase shifts (35.5) while the second 
line follows from the unitarity of S so that 5"' = S^. We can now take the limit 
/? — > CO since the R dependence has been cancelled away. 

This is essentially what we want to prove since for the left hand side we al- 
ready have the semiclassical theory for the trace of the diff"erence of Green's func- 
tions, 

d{k) - doik) = — ^-Im(tr(G(A:) - GqW) . (35.8) 
Ink 

There are a number of generalizations. This can be done in any number of 
dimensions. It is also more common to do this as a function of energy and not 
wave number k. However, as the asymptotic dynamics is free wave dynamics 
labeled by the wavenumber k, we have adapted k as the natural variable in the 
above discussion. 

Finally, we state without proof that the relation (35.7) applies even when there 
is no circular symmetry. The proof is more difficult since one cannot appeal to the 
phase shifts 5„, but must work directly with a non-diagonal S matrix. 



35.2 Quantum mechanical scattering matrix 

The results of the previous section indicate that there is a connection between the 
scattering matrix and the trace of the quantum Green's function (more formally 
between the difference of the Green's function with and without the scattering 
center.) We now show how this connection can be derived in a more rigorous 
manner. We will also work in terms of the energy E rather than the wavenumber 
k, since this is the more usual exposition. Suppose particles interact via forces of 
sufficiently short range, so that in the remote past they were in a free particle state 
labeled y6, and in the distant future they will likewise be free, in a state labeled y. 
In the Heisenberg picture the 5 -matrix is defined as S = n_Q^ in terms of the 
M0ller operators 

Q± - hm e'""^e-'"»"^ , (35.9) 

where H is the full Hamiltonian, whereas Hq is the free Hamiltonian. In the 
interaction picture the S -matrix is given by 

r— >cx> 
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Texp 



(-'£ 



dtH'{t)\ , 



(35.10) 



where H' = V = H - Hq is the interaction Hamiltonian and T is the time-ordering 
operator. In stationary scattering theory the S matrix has the following spectral 
representation 



S = 



dES{E)6{Ho-E) 



S{E) = Q^{E)Q-_\E), Q^iE) = l + iHo-E±ie)-'V , (35.11) 



such that 



Tr 



S\E)—S{E) 
dE 



Tr 



Ha 



le 



H-E-ie 



.(35.12) 



The manipulations leading to (35.12) are justified if the operators Q±(E) can be appendix K 
linked to trace-class operators. 

We can now use this result to derive the Krein-Lloyd formula which is the 
central result of this chapter. The Krein-Lloyd formula provides the connection 
between the trace of the Green's function and the poles of the scattering matrix, 
implicit in all of the trace formulas for open quantum systems which will be pre- 
sented in the subsequent chapters. 



35.3 Krein-Friedel-Lloyd formula 



The link between quantum mechanics and semiclassics for scattering problems is 
provided by the semiclassical limit of the Krein-Friedel-Lloyd sum for the spectral 
density which we now derive. This derivation builds on the results of the last 
section and extends the discussion of the opening section. 

In chapter 33 we linked the spectral density (see (31.16)) of abounded system 



d{E) = 2 d{E„ - E) 



(35.13) 



via the identity 



5{En-E) = 



1 1 
lim - Im 

e^O TT E - En + 16 

lim - lm{E„\- — — ^l^,,) 
e^O 7T E - H + ie 



^— lim {En 

2ni e^o 



1 



E-H-ie E-H + ie 



(35.14) 
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to the trace of the Green's function (34.1.1). Furthermore, in the semiclassical 
approximation, the trace of the Green's function is given by the Gutzwiller trace 
formula (34. 1 1) in terms of a smooth Weyl term and an oscillating contribution of 
periodic orbits. 

Therefore, the task of constructing the semiclassics of a scattering system is 
completed, if we can find a connection between the spectral density d{E) and the 
scattering matrix S . We will see that (35. 12) provides the clue. Note that the right 
hand side of (35. 12) has nearly the structure of (35. 14) when the latter is inserted 
into (35.13). The principal difference between these two types of equations is that 
the S matrix refers to outgoing scattering wave functions which ai^e not normal- 
izable and which have a continuous spectrum, whereas the spectral density d{E) 
refers to a bound system with normalizable wave functions with a discrete spec- 
trum. Furthermore, the bound system is characterized by a hennitian operator, 
the Hamiltonian H, whereas the scattering system is characterized by a unitary 
operator, the S -matrix. How can we reconcile these completely different classes 
of wave functions, operators and spectra? The trick is to put our scattering system 
into a finite box as in the opening section. We choose a spherical conatiner with 
radius R and with its center at the center of our finite scattering system. Our scat- 
tering potential V{r) will be unaltered within the box, whereas at the box walls we 
will choose an infinitely high potential, with the Dirichlet boundary conditions at 
the outside of the box: 



In this way, for any finite value of the radius R of the box, we have mapped our 
scattering system into a bound system with a spectral density d{E;R) over dis- 
crete eigenenergies E„(R). It is therefore important that our scattering potential 
was chosen to be short-ranged to start with. (Which explains why the Coulomb 
potential requires special care.) The hope is that in the limit /? ^ oo we will 
recover the scattering system. But some care is required in implementing this. 
The smooth Weyl term d{E; R) belonging to our box with the enclosed potential V 
diverges for a spherical 2-dimensional box of radius R quadratically, as kR^/{4k) 
or as R^ in the 3-dimensional case. This problem can easily be cured if the spec- 
tral density of an empty reference box of the same size (radius R) is subtracted 
(see figure 35.2). Then all the divergences linked to the increasing radius R in 
the limit 7? — > oo drop out of the difference. Furthermore, in the limit 7? — > oo 
the energy-eigenfunctions of the box are only normalizable as a delta distribution, 
similarly to a plane wave. So we seem to recover a continous spectrum. Still the 
problem remains that the wave functions do not discriminate between incoming 
and outgoing waves, whereas this symmetry, namely the hermiticity, is broken in 
the scattering problem. The last problem can be tackled if we replace the spec- 
tral density over discrete delta distributions by a smoothed spectral density with a 
small finite imaginary part rj in the energy E: 



m\r=R = . 



(35.15) 



n ^ 



E - En{R) -ir] E- En{R) + ir] 



1 



1 



} 



(35.16) 



scattering - 29dec2004 



ChaosBook.org versionI4.4.I, Apr 21 2013 



CHAPTER 35. QUANTUM SCATTERING 



678 



Note that d{E + iij; R) + d{E - ir]; R) = -d{E + ii]; R). By the introduction of the 
positive finite imaginary part 77 the time -dependent behavior of the wave function 
has effectively been altered from an oscillating one to a decaying one and the 
hermiticity of the Hamiltonian is removed. Finally the limit 77 — > can be carried 
out, respecting the order of the limiting procedures. First, the limit R ^ 00 has 
to be performed for a finite value of rj, only then the limit 77 ^ is allowed. In 
practice, one can try to work with a finite value of R, but then it will turn out (see 
below) that the scattering system is only recovered if /? a/?? » 1 . 

Let us summarize the relation between the smoothed spectral densities d{E + 
irj; R) of the boxed potential and d'-^^E + irj; R) of the empty reference system and 
the S matrix of the corresponding scattering system: 



This is the Krein-Friedel-Lloyd formula. It replaces the scattering problem by 
the difference of two bounded reference billiards of the same radius R which fi- 
nally will be taken to infinity. The first billiard contains the scattering region or 
potentials, whereas the other does not (see figure 35.2). Here d{E + ir];R) and 
d^^\E + iri;R) are the smoothed spectral densities in the presence or in the ab- 
sence of the scatterers, respectively. In the semiclassical approximation, they are 
replaced by a Weyl term (34. 10) and an oscillating sum over periodic orbits. As in 
(34.2), the trace formula (35.17) can be integrated to give a relation between the 
smoothed staircase functions and the determinant of the S -matrix: 



Furthermore, in both versions of the Krein-Friedel-Lloyd formulas the energy ar- 
gument E + ir] can be replaced by the wavenumber argument k + irj' . These expres- 
sions only make sense for wavenumbers on or above the real fc-axis. In particular, 
if k is chosen to be real, 77' must be greater than zero. Otherwise, the exact left 
hand sides (35.18) and (35.17) would give discontinuous staircase or even delta 
function sums, respectively, whereas the right hand sides are continuous to start 
with, since they can be expressed by continuous phase shifts. Thus the order of 
the two limits in (35.18) and (35.17) is essential. 

The necessity of the -i-/?7 prescription can also be understood by purely phe- 
nomenological considerations in the semiclassical approximation: Without the irj 
term there is no reason why one should be able to neglect spurious periodic or- 
bits which are there solely because of the introduction of the confining boundary. 
The subtraction of the second (empty) reference system removes those spurious 
periodic orbits which never encounter the scattering region - in addition to the re- 
moval of the divergent Weyl term contributions in the limit R ^ 00. The periodic 
orbits that encounter both the scattering region and the external wall would still 





-^lndet5(£'). (35.17) 




(35.18) 
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survive the first limit /? ^ oo, if they were not exponentially suppressed by the 
+iT] term because of their 



behavior. As the length L{R) of a spurious periodic orbit grows linearly with the 
radius R. The bound Rij' » 1 is an essential precondition on the suppression of 
the unwanted spurious contributions of the container if the Krein-Friedel-Lloyd 



Finally, the semiclassical approximation can also help us in the interpretation 
of the Weyl term contributions for scattering problems. In scattering problems the 
Weyl term appears with a negative sign. The reason is the subtraction of the empty 
container from the container with the potential. If the potential is a dispersing bil- 
liard system (or a finite collection of dispersing billiards), we expect an excluded 
volume (or the sum of excluded volumes) relative to the empty container. In other 
words, the Weyl term contribution of the empty container is larger than of the 
filled one and therefore a negative net contribution is left over. Second, if the 
scattering potential is a collection of a finite number of non-overlapping scatter- 
ing regions, the Krein-Friedel-Lloyd formulas show that the corresponding Weyl 
contributions are completely independent of the position of the single scatterers, 
as long as these do not overlap. 



35.4 Wigner time delay 

The term ^IndetS' in the density formula (35.17) is dimensionally time. This 
suggests another, physically important interpretation of such formulas for scatter- 
ing systems, the Wigner delay, defined as 



and can be shown to equal the total delay of a wave packet in a scattering system. 
We now review this fact. 

A related quantity is the total scattering phase shift &(k) defined as 



formulas (35.17) and (35.18) are evaluated at a finite value of R. 



exercise 35.1 



d{k) = ^Argdet(S(^)) 
dk 

= -/^logdet(S(A:) 




(35.19) 



detS{k) - e 



SO that d{k) ^ 4;@{k). 
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The time delay may be both positive and negative, reflecting attractive re- 
spectively repulsive features of the scattering system. To elucidate the connection 
between the scattering determinant and the time delay we study a plane wave: 

The phase of a wave packet will have the form: 

(f) = k ■ X - cot +0. 

Here the term in the parenthesis refers to the phase shift that will occur if scattering 
is present. The center of the wave packet will be determined by the principle of 
stationary phase: 

- dcp = dk ■ X - do) t + d& . 

Hence the packet is located at 

^ dio d@ 
x^ —t -. 

dk dk 

The first term is just the group velocity times the given time t. Thus the packet is 
retarded by a length given by the derivative of the phase shift with respect to the 
wave vector k. The amval of the wave packet at the position x will therefore be 
delayed. This time delay can similai^ly be found as 



To show this we introduce the slowness of the phase s = k/oj for which 5 • = 1, 
where Vg is the group velocity to get 

dk ■ X - s- xdoj = — doj , 



since we may assume x is parallel to the group velocity (consistent with the 
above). Hence the arrival time becomes 



doj 



V. 



'8 




+ 



da> 



If the scattering matrix is not diagonal, one interprets 
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as the delay in the j'th scattering channel after an injection in the ith. The proba- 
bility for appearing in channel j goes as |5,jp and therefore the average delay for 
the incoming states in channel / is 



where we have used the derivative, d/doj, of the unitarity relation S • S ' = 1 vahd 
for real frequencies. This discussion can in particular be made for wave packets 
related to partial waves and superpositions of these like an incoming plane wave 
coiTcsponding to free motion. The total Wigner delay therefore corresponds to the 
sum over all channel delays (35.19). 

Commentary 

Remark 35.1 Krein-Friedel-Lloyd formula. The third volume of Thirring [35.1], 
sections 3.6.14 (Levison Theorem) and 3.6.15 (the proof), or P. Scherer's thesis [35.15] 
(appendix) discusses the Levison Theorem. 

It helps to start with a toy example or simplified example instead of the general the- 
orem, namely for the radially symmetric potential in a symmetric cavity. Have a look at 
the book of K. Huang, chapter 10 (on the "second virial coefficient"), or Beth and Uhlen- 
beck [35.5], or Friedel [35.7]. These results for the correction to the density of states are 
particular cases of the Krein formula [35.3]. The ICrein-Friedel-Lloyd formula (35.17) 
was derived in refs. [35.3, 35.7, 35.8, 35.9], see also refs. [35.11, 35.14, 35.15, 35.17, 
35.18]. The original papers are by Krein and Birman [35.3, 35.4] but beware, they are 
mathematicans. 

Also, have a look at pages 15-18 of Wirzba's talk on the Casimir effect [35.16]. Page 
16 discusses the Beth-Uhlenbeck formula [35.5], the predecessor of the more general 
Krein formula for spherical cases. 

Remark 35.2 Weyl term for empty container. For a discussion of why the Weyl term 
contabution of the empty container is larger than of the filled one and therefore a negative 
net contribution is left over, see ref. [35.15]. 

Remark 35.3 Wigner time delay. Wigner time delay and the Wigner-Smith time 
delay matrix, are powerful concepts for a statistical description of scattering. The diagonal 
elements Qaa of the lifetime matrix Q - -iS^^dS/dco, where S is the [2A^x2A^] scattering 
matrix, are interpreted in terms of the time spent in the scattering region by a wave packet 
incident in one channel. As shown by Smith [35.26], they are the sum over all ouput 
channels (both in reflection and transmission) of Atab = Re [{-i/ S ab){dS at/ da))] weighted 
by the probability of emerging from that channel. The sum of the Qa„ over all 2N channels 
is the Wigner time delay tw = Qaa, which is the trace of the lifetime matrix and is 
proportional to the density of states. 



(A?,) ^ Y,\s^/AUj^Re(-iJ]s*j^) = Rei-iS'.^\ 



'ii 



J J 
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Exercises 



35.1. Spurious orbits under the Krein-Friedel-Lloyd con- 
truction. Draw examples for the three types of period 
orbits under the Krein-Friedel-Lloyd construction: (a) 
the genuine periodic orbits of the scattering region, (b) 
spurious periodic orbits which can be removed by the 
subtraction of the reference system, (c) spurious peri- 
odic orbits which cannot be removed by this subtraction. 
What is the role of the double limit 77 — > 0, container size 
b oo7 





R 



35.2. The one-disk scattering wave function. 

one-disk scattering wave function. 



Derive the 



35.3. Quantum two-disk scattering. 

classical spectral determinant 



Z(e) = Y] 



pj.i 



1 - 



A' 



j+21 



(Andreas Wirzba) 



Compute the quasi- 



for the two disk problem. Use the geometry 



The full quantum mechanical version of this problem 
can be solved by finding the zeros in k for the deter- 
minant of the matrix 



i-ir J„,{ka) 
2 Hl!\ka) 



where J„ is the nth Bessel function and H^P is the Han- 
kel function of the first kind. Find the zeros of the de- 
terminant closest to the origin by solving det M{k) - 0. 
(Hints: note the structure M = / -1- A to approximate the 
determinant; or read Chaos 2, 79 (1992)) 

35.4. Pinball topological index. Upgrade your pinball sim- 
ulator so that it computes the topological index for each 
orbit it finds. 
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